#include <iostream>

#include "sf.01.h"

int main(int argc, char **argv) {
	
	ScaleFactor a1, a2, a22;
	
	a22.z(2);
	a1.z(1);
	a2.z(0);
	
	Sec s1=0.5, s2=0.7;
	
	std::cout<< a1.H(.25, .75) << std::endl;
	std::cout<< a1.s(.25, .75, .73).Gyr() << std::endl;

	double timeUnit = 3.085678e21/1e5/1e9/365.25/86400.0;
	
	
	std::cout<< (a22.s(.25, .75, .73)).Gyr() << " Gyr"<< std::endl;
	std::cout<< (a1.s(.25, .75, .73)).Gyr() << " Gyr"<< std::endl;
	std::cout<< (a2.s(.25, .75, .73)).Gyr() << " Gyr"<< std::endl;
	
	/*
	std::cout<< (a2.s(.25, .75, .73) - a1.s(.25, .75, .73)).Gyr() << " Gyr"<< std::endl;
	
	std::cout<< (a2.s(.25, .75, .73) - a1.s(.25, .75, .73)).Gyr()/timeUnit/300 << " MaxTime"<< std::endl;
	*/
	
	
	
	
	
	return 0;
}

